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Abstract. In this paper we investigate the behavior of a family of steady state solutions of a nonlinear 
reaction diffusion equation when some reaction and potential terms are concentrated in a e-ncighborhood of 
a portion T of the boundary. We assume that this e-neighborhood shrinks to T as the small parameter e goes 
to zero. Also, we suppose the upper boundary of this e-strip presents a highly oscillatory behavior. Our main 
goal here is to show that this family of solutions converges to the solutions of a limit problem, a nonlinear 
elliptic equation that captures the oscillatory behavior. Indeed, the reaction term and concentrating potential 
are transformed into a flux condition and a potential on T, which depends on the oscillating neighborhood. 



1. Introduction 

In this work we analyze the behavior of a family of steady state solutions of a homogeneous Neumann 
problem for a nonlinear reaction diffusion equation when some reaction and potential terms are concentrated 
in a e-neighborhood of a subset T of the boundary that shrinks to T as the small parameter e goes to 
zero. Roughly, we are assuming that some reactions of the system occur only in an extremely thin region 
near the border with oscillating upper boundary. We show that in some sense this singular problem can 
be approximated by a nonlinear elliptic system with nonlinear boundary conditions where the oscillatory 
behavior of the neighborhood is captured as a flux condition and a potential term on the portion T of the 
boundary. 

To describe the problem, let £1 = (0, 1) x (0, 1) be the open square in R 2 and let T c dfl be the line 
segment in R 2 given by 

r={(i,0)£R 2 : x 6(0,1)}. 
We consider the following uniformly bounded e-neighborhood of T with oscillatory upper boundary 

uj € = {(x,y) £ R 2 : x £ (0,1) and < y < e G e (x)}. 

Here we assume that G e (-) is a function satisfying < Go < G £ (-) < Gi for fixed positive constants Go and 
Gi which oscillates as the small parameter e — > 0. This is expressed by assuming that 

G e (x) = G(x,x/e). (1.1) 

The function G : (0, 1) x R i— > R is a positive smooth function, with y — > G(x,y) periodic in y for fixed x 
with period l(x) uniformly bounded in (0, 1), that is, < ^o < < Let us observe that our assumptions 
includes the case where the function G £ presents a purely periodic behavior, for instance, G e (x) = 2+cos(x/e). 
But it also considers the case where the function G e defines a strip where the oscillations period, the amplitude 
and the profile vary with respect to x £ (0, 1). See Figure [l] below that illustrates the oscillating strip w c C il. 
We are interested in the behavior of the solutions of the nonlinear elliptic problem 

-Au e + \u e + -X Uc V € (-) u e = -X Ue /„(-, u e ) + /i(-, u e ) in Q 

= onffl 

on 
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Figure 1. The open square ft and the e-strip w e 



where X Uc is the characteristic function of the set u) e , n denotes the unit outward normal vector to dfl and 
A > is a suitable real number. The nonlincaritics /o, /i :0xli->K are continuous in both variables and 
C 2 in the second one, where O C M 2 is an open set containing f2. We also assume that there exists C > 
independent of e such that the family of potential V e G L°°(Q) satisfies 

1 



\V e (x,y)\ z dxdy<C. 

Also, we suppose there exists a function Vq G L 2 (Y) which is the limit of the concentrating term 

1 



lim 

e-S-0 e 



V (pdS, \Jlp eC°°{VL). 



(1.3) 



(1.4) 



In our model, we use the characteristic function X L 
concentration on the region u e C by the term 

1 



and the small positive parameter e to express the 



x Uc e L°°(fi). 

Since u> e C (0, 1) X (0, e G\) is thin and it is "approaching" the line segment T C 9f2, it is reasonable to expect 
that the family of solutions u e will converge to a solution of an equation of the same type with nonlinear 



boundary condition on V. Indeed, we will show that, for A big enough, the solutions of (1.2) converge in 
to the solutions of the nonlinear elliptic problem 

—Au + Xu = fi(-,u) in f] 

du 



dn 



+ V (-)u = n(-)f (;u) onT 



(1.5) 



du 
dn 



= on an \ r 



(1.6) 



where the boundary coefficient /i £ L°°(r) is related to the oscillating function G e and is given by 

= Ti~\ / G ( x ,y)dy, Vac e (0,1). 

As previously mentioned, we obtain a limit problem with nonlinear boundary condition that captures the 
oscillatory behavior of the upper boundary of the set tu e . In fact, its nonlinear boundary condition involves 
the function fi(x), that is, the mean value of G(x, •) for each x € (0, 1). We now summarize these assertions 
precisely as our main result. 



Theorem 1.1. Let u e be a family of solutions of the problem (1.2) satisfying \\u e \\ L oo^ < R for some 
positive constant R independent of e. Then, there exists A* G K independent of e such that, for all A > A* . 
there exists a subsequence, still defined by u e , and a function u G H 1 ^), with ||it||j>o(n) < R, solution of 
the problem (1.5) satisfying \\u e — u||jji(n) — > as e — > 0. 



TERMS CONCENTRATING ON BOUNDARY 



3 



Since we are concerned with solutions which are uniformly bounded in L°°(Q), we may assume that the 
nonlinearities /o and f\ are C 2 -functions with bounded derivatives in the second variable. Indeed, we may 
perform a cut-off in /o and fi outside the region \u\ < R, without modifying any of these solutions. 

This kind of problem was initially studied in [4] where linear elliptic equations posed on C 2 -regular domains 
were considered. There, the e-neighborhood is a strip of width e with base in a portion of the boundary 
without oscillatory behavior. Later, the asymptotic behavior of the attractors of a parabolic problem were 
analyzed in [5J[3], where the upper semicontinuity of attractors at e = was proved. The same technique 
of [1] has been used in [TJ 12] , where the results of [H [5] were extended to reaction-diffusion problems with 
delay. 

The goal of our work is to extend some results of [4] to nonlinear elliptic problems when the upper boundary 
of the strip ui e presents a highly oscillatory behavior. We also use some ideas of [21 [51 [5] where elliptic and 
parabolic problems defined in thin domains with a highly oscillatory boundary have been extensively studied. 
It is important to note that in our work the boundary dfl of the domain f2 is only Lipschitz . 

2. Some technical results 

In this section we describe some technical results that will be needed in the proof of the main result. We 
initially analyze how our concentrating integrals converge to boundary integrals. We adapt the results of [1] 
on convergence of concentrated integrals using that < Go < G e (-) < G\ uniformly in e > 0. 

Lemma 2.1. Suppose that v G H s (n) with 1/2 < s < 1 and s — 1 > —1/q. Then, for small e , there exists 
a constant C > independent of e and v such that for any < e < €q, we have 



\ I i«i a de<c||tc. 



(O) 



Proof. Note that 

1 /",.„„ 1 



/ \v\ q dt< - f \v\ q d^ 



e 

where r e is the strip of width e G\ and base in Y without oscillatory behavior given by 

r t = {0,2/) e M 2 : x e (0,1) and < y < ed} . 
Since u;, is contained in r,, the result follows from 4, Lemma 2.11. □ 



As a consequence of Lemma 2.1 and [A, Lemma 2.5] we obtain 



Lemma 2.2. Suppose that the family V £ satisfies fl.Sfy and (1-4). Then, for s > 1/2, a > 1/2 and 



s + a> 3/2, if we define the operators T e : H S (Q) -> (H a {n)Y by 

(T e (u),(p) = - I V e u(pdt;, for e > 0, and (To(u),ip) = I VoU(pdS, 



e _ 

we have T e -> T in C(H s (n), (H a {Q))'). 

We also need some results about weak limits of rapidly oscillating functions. 



r 



Lemma 2.3. If G e is defined as in (1.1), then 



1 



Ge(-) -> K-) = J o G(; S )ds W*-L°°(P,i). 



Proof. We have to prove 

K 




— / G(x, s) ds > <p(x) dx as e -t 0, Vw G ^(0, 1). 
x ) Jo J 
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Since the set of step function is dense in T/^O, 1) and any step function is a linear combination of characteristic 
functions, it is enough to take <p = X( e ,f)i the characteristic function on the interval (e, /), for < e < / < 1. 
So, we have to estimate the integral 



— r / G(x, s) ds > dx. 
x ) Jo J 



For this, let r\ > be a small number and let {e = Xq, x±, x n = /} be a partition for the interval (e, /), 
and Xi be a fixed point in the interval Ji — [xi-\, Xi], i — 1, n, such that 

sup sup \G(x,y) - G(xi,y)\ < rj. 
Observe that we can write I e j = 2j=i /) where 

n p 

T lf =J2 {G{^x/e) - G(xi,x/e)}dx, 

i=l J J i 

?e,f = / , J { G(xi,x/e) - ^-y J G(xi,s)ds > dx 



•I, - 1 1 { jf ' ■) * - of* •) *} * 



G(cc, s) ds > da; 



It is easy to estimate the integrals 7* y, 7|y, I* t and 7|, to obtain |J^.| < Tj(f — e), |7f *| < vif ~ e ): 
1^/1 < Gi ||?" - i||i-(o,i) (/ ~ e) and ^ < d (h/l 2 ) ||P' - i|U~ (0j i) (/ - e), where the function f" is 
the step function defined for each 77 > by ^(a;) = Z(iEj) as a^ G Ji. Since these inequalities do not depend 
on e > 0, and \\l v — /||_£,oo( ,i) as r\ — > uniformly in e, we have that 7* y, 7|y, 7* * and 7^ goes to zero 
as — ?> uniformly in e > 0. Hence, to conclude the proof, we just evaluate the integral 7| But this is a 
application of 7J Theorem 2.6]. □ 

The following result will also be needed. 
Lemma 2.4. Suppose that h, <p 6 H s (fl) with 1/2 < s < 1. Then, 

lim- / htpd£= [ fj,-/(h)j{cp)dS, (2.1) 



e->0 e 



r 

2/ 



where /1 € 7°°(r) zs given by (1.61 and 7 : H S (Q) t— > L (T) is i/ie trace operator. 
Proof. Initially, let ft and </? be smooth functions defined in f2 independent of e. Note that 



If. /" . .„ 1 



e ./,.,_ ./r e 



1 (.eG.ts) /•! 



*/ 



hipdt;— / fihtpdS=- / / h{x,y) tp(x,y) dydx — I /j,(x) h(x,0) ip(x,0) dx 



Taking y = eG c (a;) z, we obtain by changing variables 



1 r eG t (x) pi pi 

h(x,y) ip(x,y) dydx — / / h (x, eG e (x)z) if (x, eG e (x)z) G e (x) dzdx. 



^0 ^0 ^0 
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Thus, adding and subtracting h(x, 0) ip(x, 0) G e (x) dx, we get 
hipdt; — / fihtp dS 



< 



h(x, 0) ip(x, 0) G e {x) dx — I fj,(x) h(x 7 0) ip(x,0) dx 

Jo 

I. r l 

G e (x) [h (x, eG e (x)z) ip (x, eG e (x)z) — h{x, 0) tp(x, 0) ] dzdx 



Now, since e G e (a:)z — >■ as e — ?► uniformly for (x, z) <E [0, 1] x [0, 1], we have 



G e (x) [h {x, eG t (x)z) ip (x, eG e (x)z) — h(x, 0) <p(x, 0)] dzdx 



as e — > 0. 



(2.2) 



Hence, we obtain from Lemma 2.3 and (2.2 1 that 



hipdt; — / nh<p dS 



as e — s- 0. 



Consequently, the proof of equality (2.1 1 follows by density arguments and the continuity of the trace operator 
7 (see [TD]). □ 

3. Abstract setting and existence of solutions 



To write the problems (1.2) and (1.5 1 in an abstract form, we first define the continuous bilinear forms 
a e : H 1 ^) x if 1 (ft) H> R, e > 0, by 

a e (u, v) = [ WuWvdxdy + X j uvdxdyA — / V^uvdxdy, for e > 0, 

do (u, v) = / Vu Vv dxdy + A / uv dxdy + / Vquv dxdy. 
Jn Jn Jr 

We then consider the linear operator A e : H 1 C -ff _1 (f2) h-> i/ _1 (f2) defined by the relationship 

(A £ u,w) _ 1 1 = a e (zt, u), for all v € /f^fi). 



Now, we can write the problem (1.2) as A e u — F e (u), for e > 0, where F t : 1 (SI) h-> s (f2) with 
1/2 < s < 1 is defined by 



(F 0l e(u),«) = ^ / /„(;. and =//.,(;. »)r<. tor all r .£//-( O). 



(3.2) 



Similarly, we can write the problem ( 1.5 ) in an abstract form as Aqu = F(u), where F : i/ 1 (Sl) i-> H S (S1) 
with 1/2 < s < 1 is defined by 

F = Fq + F\ , 

Fi is given by dOk and (fb(u),t;) = / M7(/o(£, «)) 7(«) for all v e tf s (Q), 



where <E L°°(r) is given by (1.6) and 7 : H s (fl) H> L (T) is the trace operator. 

Lemma 3.1. There exists A* € K independent of e > swc/i £/iai £/ie bilinear form a t is uniformly coercive 
for all A > A* . 

Proof. Here we will just consider the case a e with e > 0. A similar argument gives the result to the bilinear 
form qq. First we note that for every G iJ 1 (S7) we have 



ae(^^)>||V0||l 2(n) + A||0||^ (n) -- / (V e )_ \4,\ 4 dt 



(3.4) 



(i 
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where (V e )_ is the negative part of the potential V e such that V e = (V 6j 
have the following bound 



(V e )_. For this negative part we 



(V e )_ \0\ 2 d£ < 



m 2 « 



\<t>\ <%) 



(v e )_\4>\ 2 d^cu\\ 



,2(1- S ) 
\L 2 (Q) 



Choosing | < s < 1 and s — 1 ^ — 4, that is, § ^ s < 1, and using the Lemma 2.1 with g = 4, we get 

1 

e 

Next we can use Young's inequality to obtain 

1 
e 



(Ve)-\ct>\ 2 d^6\ml 1(n) +C s \\c!>f 



L 2 (Q) 



(3.5) 



for any 5 > 0. Then, it follows from (3.4 1 and ( |3.5| ) that 

a £ (0, 0) ^ (A — (1 + C,)) |H|* 2(n) +(1-5) U\\ 2 HHa) . 

Consequently, we can take 5 > small enough and A > large enough such that 

a e (<j>,<f>)>C\\cf>\\ 2 HliQ)) V^eH^n), (3.6) 

with C — C(A) > independent of e. Therefore, the bilinear form a e , < e ^ eo, is strictly coercive. We 
still note that if V e ^ in (3.1 ), then we can take any A > 0. □ 



The Lemma 3.1 implies that u e is a solution of Ok if only if u € £ H 1 ^) satisfies u c = A^F^u*), that 
is, u e must be a fixed point of the nonlinear map 

A' 1 oF t : H l (Q) i-> H 1 ^), 



for all A > A*. The existence of solutions of (1.2) follows then from Schauder's Fixed Point Theorem. In a 
very similar way, the solutions of the limit problem (1.5) can be obtained as fixed points of the map 



Note that they can also be obtained as limits of solutions of (1.2), as shown in the next section. 

4. Upper semicontinuity of steady states 

In order to obtain the upper semicontinuity of the family it e , we study the behavior of the maps F t and 
F defined in OJ and fl33| . 

Lemma 4.1. (1) If u £ H 1 ^) satisfies \\u\\ L<*>(n) < R> then there exists K > independent of e such that 
sup{||F(«)|| fl _. (n) ,||F e («)|| H _. (n) } <K, for all 1/2 < s < 1. 

(2) Suppose ||it||fl-im)nz,~(r2) < R- Then, we have for all 1/2 < s < 1 

||.F e (it) — F{u)\\h-s(q) ~^ as e 0, uniformly in u. 

(3) .4Zso, ifu'^u in H l (£l), then ||F e (u £ ) - F(u)||^_ s(0) -> as e ->■ 0. 



Proof. We can get (1) from continuity of the nonlinearities and Lemma 2.1 Part (3) follows from (2) adding 
and subtracting F e (u). So, we just have to prove (2). For this, it is enough consider the maps -Fo.e an d Fq 
given by (3.2) and (3.3). Initially, take Sq satisfying 1/2 < sq < 1. For each u £ H l (Q) and <j> £ H S °(H,) 



\(F 0te (u),<f>)-(F (u),<t>)\ = 
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From Lemma 2.4 we have (F 0j£ (u), <t>) — > (F (u),(f>) as e — > 0, for each (/> G _ff s °(£l). Also, this limit is 
uniform for (f> in compact sets of H s °(£l) since {Fo ]£ (w) G H~ s "(£l) : e G (0, eo]} is equicontinuous for fixed 
it G iJ 1 (rj). Hence, for 1/2 < sq < s < 1, we have that the embedding H s (il) H S( '(£l) is compact, and 



\FqA u ) ~ F o(u)\\ H ->(n) 



sup 



\(F 0>e (u) - F {u), 0)| ^ 



as e 



0. 



(4.1) 



Now, we show that (4.1 1 is uniform for u G with IMI^i^ < R. First, observe that F^ :<Ll Fq : h-> 

H~ s (il) are continuous in with the weak topology, for e > 0. Indeed, it follows from Lemma 2.1 and 

regularity of /o that there exist C, K > 0, independents of e, and < 9(x) < 1 with x G ft such that 



\F , e (u) - F , e (v)\\ H _ sm < C 



\d u f (x,9(x)u(x) + (1 - 6{x))v{x))\ 2 \u{x) - v(x)\ 2 dx 



Similarly, we get 



< CK \\u - v\\ Hs{a) , Vu,veH s (Q) 



\\F (u) - F (v)\\ H . B(n) <CK\\u- v\\ Hs(n) , for all u, v G H S (Q). 



So, since if 1 (ft) <^-> iJ s (ft) compactly for s < 1, we have that Fo j£ and are continuous in iJ 1 (f2) with 



weak topology. Hence, F 0i£ and F are uniformly continuous in {u G 1 (f2) 
proving the result. 

4.1. Proof of the main result. 



< R} c ff 1 ^), 
□ 



Proof. Now, we are in a position to prove our main result, Theorem |1.1| First, we prove that a family of 
solutions u e of (1.2 1 satisfying ||u e || i oo( f2 ) < R is uniformly bounded in i? x (ft) with respect to e. In fact, we 
have that u e satisfies 



a e (u e ,<p) = (F , e (u e ),y>) + (^i(u £ ),(^) , V^Gtf^ft) and e > 0. 



(4.2) 



Hence, if we take if = w € , we obtain from Lemma 4.1 that |a e (u e ,u e )| < i4T||ii € ||#i(n)- Due to Lemma 3.1 
we have that the bilinear form a £ is uniformly coercive for e > (we are using A > A*). Thus, 

c|MI 2 ffl(0) < k«0| < K\\u e \\ HHn) , 

for some constant C > independent of e, which implies ||M e ||ifi(o) < K/C for all e > 0. 

Since u e is uniformly bounded in if 1 (ft), we can extract a weakly convergent subsequence, still denoted 
by u e , such that u e — * it, to — if 1 (ft), for some it G -ff^ft). It is easy to see that it satisfies our limit problem 
(1.5|. Indeed, we use Lemma 2.2 and Lemma 4.1 to pass to the limit in (4.2) as e — > 0, and obtain 

a (u,<p) = (F (u),(p) + (Fi («),¥>) , for all ip e if 1 (ft). 

Now, we can prove that it e — > u in ii^(ft) showing the convergence of the norms. For this, we pass to the 
limit in a e (u e ,i/) = (F ^(u e ),u e ) + (Fi(u e ),u e ) and we use that the norm is lower semicontinuous with 
respect to the weak convergence (see [7J Proposition 1.14]), that is, ||tt||jji(n) < liminf ||u e || ^l^) , to obtain 



Vu\ z d£, < liminf / |Vu e | 2 d£ < limsup / \S7u € \ z d£, 



,e|2 



< 



\\u\'d£- / V r |4 , d£+ / iifoti,u)udt+ / m,u)u<% 



[ \Vu\ 2 d£. 



a 
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5. Final conclusion 

We have shown that the steady state solutions of a homogeneous Neumann problem for a nonlinear 
reaction diffusion equation converge to a certain limit problem when some reaction and potential terms are 
concentrated in a small neighborhood of of the boundary. 

In our analysis we showed that the family of steady state solutions converges in ii^-norm to a solution u 
of an equation of the same type, with a nonlinear boundary condition that captures both the profile and the 
oscillatory behavior of the boundary. 

An important feature here is that we are dealing with the case where the boundary presents a highly 
oscillatory behavior and, as a consequence, the limit problem is not obvious from the start. It is also worth 
remembering that the domain considered here is only Lipschitz continuous. 

We use some results developed in [3] that describe how different concentrated integrals converge to surface 
integrals to obtain a rigorous strong convergence result in Theorem 

A natural question is whether such approximation results can be improved in order to describe the 
asymptotic behavior of the Dynamical System generated by the parabolic equation associated to the problem 



(1.2 1 posed in more general regions of the R . It is our goal to investigate this question in a forthcoming 
paper. 
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